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Abstract 

A mathematical model of radiotherapy is proposed. The study used the classical 
24 hours way of fractionation with a weekend pause. We introduce the matrices 
of "radiotherapy" and "growth". We developed an equation of the fraction cell 
evolution, which we solved numerically. The results indicate that the accelerated 
growth of cells occurs due to the decrease of the fraction of slowly growing cells and 
increase of the cells that are fast growing. 



Repopulation is a serious problem in cancer radiotherapy. The growth of tumor in 
final stage of radiation can prevent the healing of a patient. Many studies have shown 
the importance of timing in radiotherapy [2]. If there is a delay in the treatment caused 
e.g. by interrupts in radiation then this the increases time given to cells for accelerated 
growth. It was shown that breaks in radiation especially after fourth weeks of treatment 
lead to the worse results, approximately 2-4.8 percent growth per day of delay. 

Why such a strange phenomenon occurs? Tumor fights for its live and it behaves 
according to the Lenz rule: the lower number of cells the faster the growth. According 
to Trott [3] accelerated repopulation occurs when the number of cells in tumor decreases 
below 1000 cells. It can be explained in a few different ways. 

Referring to a recent study there [4] are three theories for explaining reasons for 
repopulations: 

1. Fowler model, in which the author claims that cancer volume doubling time Td 
approaches potential time T pot as a result of loosing cells. 

2. Jones Model which proposes the following explanation: the tumor possesses sub- 
populations of cells growing with different velocities (speeds); cells are dying equally but 
those dividing faster gain advantage during breaks in radiation. 

3. Trott-Kummermehr model which can be called Dragon Theory. Like a mediewal 
knight cutting dragon's head have met next two new heads, here stem cells switch from 
asymmetrical division to symmetrical one. At each division from one stem cell two are 
arising. 
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1 Assumptions 



We assume as a basis the Jones model. Cancer tumor is heterogenic; it means that there 
exists fractions of cells that differ in terms of access to oxygen or nutrition or number of 
mutations. We assume there are three fractions of cells in the tumor, which we will denote: 

Xo - small number of mutations and low growth velocity 

X\ - intermediate number of mutations and medium growth velocity 

x 2 - large number of mutations and fastest growth velocity 

where these variables are normalized by 

2 

I>i = 1 

i=0 

In individual fractions there is a well determined number of cell, where 
Hi number of cells belonging to the i-th fraction 

1=1 

and where N is a total number of cells. 

One of mutation factors is the radiation itself [1]. There is no reason to prevent 
such a phenomenon during radiation. Subsequently to the radiation of tumor after each 
consecutive dose, the number of cells in tumor will decrease and cells in each individual 
fraction will undergo mutations. Below is a new model of decreasing of number of cells 
in each fraction. 

Fraction x° can be expressed by 

*°i = Vi/N (1) 

Growth of tumor is a result of growth of individual fractions of tumor cells. Each 
fraction x^, x®, x® grows with its own velocity v , v±, v 2 . Time of duplication of tumor 
T dyi determines the velocity according to: 

^ = -^- L ( 2 ) 

J-d,i 

where is volume doubling time for individual fraction tumor. There is besides velocity 
the influence on the tumor growth, on the number of cells in the particular fractions of the 
tumor. Average velocity of the tumor growth can be described by means of the formula: 

2 

$ = ViXi (3) 

i=0 
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2 Matrix of radiation 



The decrease of the tumor volume, i.e. 
quadratic formula: 



waste of stem cells is described by the linear 



N = N e~ ad - pd2 



(4) 



where N is the number of survive cells that radiotherapy , and N is the initial number 
of cells. Another form of this formula is following: 



S = e 



-ad-/3d 2 



(5) 



where S is fraction of surviving cells, a, f3 are coefficients, d is a radiation dose. 
Accordingly 

S = N/Nq (6) 



or 



s = E*J i = E*o (7) 

i i 

Here x° are initial fractions before radiotherapy and x\ after first dose of radiation. 

D^D^ (8) 



When move from fractions to the number of cells then the equations take the form: 



(9) 
(10) 

(11) 



In matrix notation we have: 



-j— ad— /3d 2 




Vo 
Hi 

\yl 



(12) 

WiJ V u ° 

The matrix is diagonal and it shows that each fraction decreases according to the 
linear-quadratic formula and there is no exchange of cells between individual fractions. 

This description suggests that all cells behave the same way and are equally sensitive 
to the absorbed dose gained by the tumor. Investigations show that cancer tumor does not 
possess uniform cells the individual cells differ in access to the oxygen or nutritious means. 
One of mutation factors is the ion radiation. During radiation surviving cells inherit 
improved conditions of oxygenations and nutrition and are undergoing rapid mutation. 
All these changes lead both to the decrease in the number of cells in individual fractions 
and also to the change of the proportions of individual fractions. We introduce coefficients 
Q and P to describe the probability that the cells from fractions Xq, x\ will shift to fraction 
x\ and X®, respectively. 
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Equation describing this process have the following form: 

vh = y°o(e- ad - fitp - Q) (13) 
y\ = y°de- ad -" d2 - P) + Qy° (u) 

y\ = y° 2 e- ad -^ + Py« (15) 

(16) 

Situation after n steps is described by the following equations: 

yt ] = (e^ d +^-Q)) yt 1] (17) 

y{ n) = Qyt^ + ie-^^-Pijy^ (18) 

y in) = P{y ( r ^ +y (n-l) e - {ad+m (1Q) 

(20) 

or in matrix notation: 




n 


fy°o\ 








\yl) 



e -ad-pd* _ Q 

Q e ~*d-f3d* _ p o | | | (2i) 

Q p e -ad-/3d 2 

We will call the matrix appearing above a radiation matrix and we will denote it 71: 

/ e -ad-P#_Q o V 

71= Q e -*d~f3d* _ p o | ^2) 

V P e -ad-f3d* J 

It can be shown by induction that for each n after summing up rows we obtain following 
equations: 

N = N e ~ n{ad+pd2) (23) 

It means that the matrix 7Z describes the diminishment of tumor cells after radiation 
according to the linear quadratic form. Additionally it shows how individual fractions 
change in time during radiation. Coefficients Q and P allows exchange of cells between 
individual fractions. 

Inserting here values for a, (3, P, Q the number of cells N and values of fractions 
Xq, X\ we can calculate the rate of decrease of the number of cells in each fraction. From 
computer simulations it follows that the vector (xi,x 2 ,x 3 ) tends to the equilibrium state 
(0,0, x 2 ). From this we conclude, that radiation of the tumor leads to the selection of 
cells which are the most mutated and which grow with the largest speed. In the limit of 
large n this equation has the form: 

$ = v 2 x { 2 n) 

It means that speed of tumor growth is larger when the tumor diminishes, and this 
effect cannot be avoided. 
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3 The growth matrix 



We describe the rate of tumor growth according to the Sole [5]. We make use of the 
results of the paper [5] in which it was shown that the equations for growth of the cell 
fractions can be written as: 



^7T = v x (l-Q') - x &(xo,Xi,x 2 ) 
at 

dx\ 

— = ViXi(l - P') + VqXqQ' - Xi$(x ,Xi,X 2 ) 

at 



dx 2 
~dt 

or in the matrix notation: 



= V 2 X 2 + ViXiP' - X2®(xo,Xi,X 2 ) 



X 



M~x* 



(24) 



where af is the time derivative of the of 
mixing matrix M. is given by 



M = 



( foil -Qf) - ${x ,x u x 2 ) 

foQ' 
V o 



Xi : x 2: x s ) , 1 is the identity matrix and 









f 1 (l-P')-$(x ,x 1 ,x 2 ) 

frP' f 2 -<5>{x ,x u x 2 ) 



(25) 

and we denote this matrix M. as the growth matrix. Equations of Sole describe how 
ratios of fractions are changing during radiotherapy. However growth of the cell number 
in each fraction which occurs during the pauses in radiotherapy we obtain through the 
following procedure: The result of the Sole equation expressed in fractions Xi we convert 
to the integer valued number of cells yi (see below). Obtained number of cells we multiply 
by factor e ln( - 2 ^* Vi = 2 Vi and next we pass from the number of cells back to fractions: 



fy'o 

\V2. 










':) 







(26) 



Here y[ is the number of cells after division during pause between consecutive pulses of 
radiotherapy. Let 

/2 V0 \ 



V = 





V o 



2 vi 







2 V2 



(27) 



denote the division matrix. The growth matrix we construct from mixing and division in 
the following way: 

g{y i ,x i )=V(y i )M(x i ) (28) 
We introduce the dependence of v 2 on v\. 



v 2 = a Vl TP(9,Q,P,Q ',P', d,n). 



(29) 



5 



where 

iP(e, Q, P, Q', P', d, n) = e e^VQ^*+y/Q*+r** (30) 
This form is for radiation period, during weekend the form is different: 

^(0, Q', P', d, n) = e e-n(VQ' 2 +?' 2 (31) 

It contains a threshold after crossing this threshold velocity v 2 decreases, in accordance 
velocity $ also diminishes. In the paper by Sole at al [5] it is shown for which parameters 
Q', P', M' the fraction x 2 exists. 



4 The radiation 

The act of radiation consists of: 

a. time of radiation — pulse radiation 

b. time between consecutive exposures 

These equations describe growth of tumor from the moment of ending of radiation till 
the next exposures. In classical fraction this time is 24 hours. The radiation lasts for a 
very short period of time: a few minutes, in comparison to 24 hour waiting period between 
consecutive fractions. Exposition after Wheldon [6] we can call pulse radiotherapy. In 
the remaining time, tumor cells repair damage from radiotherapy and they divide. The 
tumor growth appears. 

In classical radiotherapy we radiate once a day during the 4-7 weeks, depending on the 
radiation dose quantity. Duration of radiation is very short (a few minutes), remaining 
time is spent on the repair of post-radiation damage. Symbolically we can demonstrate 
this in the following form: 



(growth radiation) 5 growth 2 . . . growth 2 (growth radiation) 5 \N >= 
((growth radiation) 5 )(growth 2 (growth radiation) 5 )™^ 1 \N > 

where n is the number of weeks of radiotherapy and No is initial number of tumor cells, 
i.e. \N >= \initial number of tumor cells > 
Radiation is represented by the matrix: 

, e -ad-(3d 2 _Q o 

TZ = Q e -ad-^ _ p o | (32) 

V P e 

and 




M = 



( f (l-Q')-$(x ,x 1 ,x 2 ) 

foQ' 
V o 



h(l-P> 





- §(x ,x 1 ,x 2 ) 

hp' 








f 2 - $(x ,x 1: x 2 ) J 



(33) 
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The radiation process can be described by the following equation: 



y (n+l) = Uy n ( 34 ) 

This equation describes the diminishing of the number of cells in fractions after a pulse 
of radiation. Next we can write the equation which describe the rise of tumor until 
next exposition. First we describe the change in proportion in fractions in tumor cells. 
To this aim we change the variable from y to x. We assume that the proportion of 
y(t + next day)/y(i) is the same as normalized variables: 

y(t + next day) x(t + next day) 

W) = W) { } 

The Sole equations describe change of proportion during growth of tumor between expo- 
sitions: 

^— = Mx n (36) 
dt v 1 

Differential equations were solved on the interval of one working day and on the interval 

of three days during weekend. We change solutions of this equation again into the number 

of cells. Then we calculate how many new cells arise during the pause between radiations. 

We again pass from x back to y as before and we use following equation: 

y> n+1 = VT l+1 (37) 

We repeat this procedure during prescribed cure time expressed in weeks. 



5 Results 

We performed numerical calculations for two different sets of parameters: in the first 
case the probability coefficients of Q, P, Q', P' were zero, in the second they were different 
from zero and were Q = 0.0005, P = 0.0005, Q' - 0.1, P' = 0.1. We took the parameters 
a = 0.2,(3 = 0.02, d = 2Gy,n = 30, v = 0.01, Vi = 0.016. v 2 was calculated from the 
formula ?? for m = 5 and threshold 6 = 0.005. At Q, P, Q', P' equal zero velocity of 
growth rising of tumor was initially slower, and next faster. The conclusion is simply: 
the fact of existence of the population of cells of with different growth rates is sufficient 
for the tumor to grow faster as the corollary of the radiation and intervals between the 
fraction radiotherapy. At this point the faster growing cells are gaining the population 
dominance over the slower growing cells and finally lead to the accelerated growth tumor. 
After introduction the P, Q, Q', P' different from zero the velocity of the tumor growth 
was also larger. According to this model probabilities Q, P, Q', P' are responsible for 
the velocity change via the change of the cell population x from x Q and X\ to x 2 - Fraction 
x 2 is the fastest growing. Biologically it can be explained that fraction x 2 gains the 
best condition for growth due to the improvement of oxygenation and better nutrition. 
Additionally fraction x 2 is built from the most undifferentiated cells and most mutated 
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cells. The factor that causes the decrease fractions xq and x\ and increase of the fraction 
X2 is the radiotherapy which leads to the death of cells and simultaneously by causing 
the shift of mutation from xq and X\ to x^. In this way additional radiation influence 
leads to some differences between situation when the ionization energy is the cause of 
destroying only cells in fractions (Q — 0, P — 0, Q' — 0, P' = 0) and leads to the shift of 
cells from one fraction to another (Q ^ 0, P ^ 0, Q' ^ 0,P' ^ 0). These differences can 
be seen after subtraction of two graph, see Figs. 1 and 2. In this way we can see how 
the change of the oxygenation, better nutrition and mutation influence growth velocity. 
According to the assumption about the existence of the threshold the influence of the 
coefficients (Q ^ 0, P ^ 0, Q' ^ 0, P' ^ 0)is as following: In the first stage we see that 
they accelerate tumor growth and next they cause the slowing down of the tumor growth. 
On this example we see a new mechanism leading to the death of the tumor. Namely 
the existence of the 9 threshold, Namely the existence of the theta threshold, causes that 
growth mutacji not only shifts cells to faster growing fraction but also crossing sufficient 
threshold causes slower growth of tumor. 

6 Conclusions 

Our model explains accelerated growth tumor according to Jones model. The existence 
of cells fraction of different growth velocity is the cause of the accelerated growth tumor. 
During the radiation therapy because of the interrupts, cells have time to take advantage 
of the differences in the velocities of growth and increase the number of cells in fastest 
growing fraction. In the course of radiation the living conditions of cells are changing: the 
oxygenation and nutrition is better and additionally mutations appear. Together these 
factors causes the change of the number of cells in the particular fractions. We have 
described these changes by the coefficients P and Q responsible for the changes during 
the radiation pulse and coefficients Q', P' responsible for the changes occurring between 
different fractions. These coefficients are modifying the tumor growth: in the begin- 
ning they accelerate and after crossing the threshold O they slow down. The existence 
of the threshold O could explain the benefits from the simultaneous radiochemotherapy. 
Chemotherapy has the mutagen function. This mechanism in conjunction with radiother- 
apy facilitates the crossing of the threshold O after which the tumor growth is slowing 
down. From the model it follows that natural state of the tumor is the state described 
by the vector (0,0,1), it means that it tends to the fastest fraction. The fastest and 
least differentiated fraction gains the crucial dominance. It seems to be in accordance 
with the clinical experience that often the revival of the tumor is more malicious and less 
differentiated. 

Model that we introduced is the enlargement of the linear-quadratic model. When 
we take into account only total number of cells it describes the diminishing of the cells 
tumor exactly the same way as the linear-quadratic model. However it allows to see how 
the numbers of cells in each fractions change during the radiation and which influence on 
the velocity of the tumor growth is the appearance of the threshold O mutation. 
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It is possible to apply our model to arbitrary doses, time of radiation and breaks in 
radiations time. 
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Figure 1: Plot of velocity of the tumor growth as a function of time (days) of radiotherapy 
Here Q,P,Q',P' are zero. 
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Figure 2: Plot of velocity of the tumor growth as a function of time (days) of radiotherapy. 
Here Q, P, Q', P' are different from zero. ^ 
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Figure 3: Plot of velocity for zero valued parameters subtracted from velocity for param- 
eters different from zero. 
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